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Abstract. We present a new interpretation of the ultradiscretization procedure in 
terms of a non-archimedean valuation on the field of Puiseux series. We shall demon- 
strate how the new interpretation of the ultradiscretization gives a correspondence 
between singularity confinement for discrete mappings and ultradiscrete singularity 
confinement for ultradiscrete mappings using results from tropical geometry. 



Integrable cellular automata are dynamical systems that are discrete in time, space and 
state that possess properties typically associated with integr ability [17]. These properties 
include the existence of Lax-pairs, the existence of conserved quantities and key properties 
regarding singularities such as the Painleve property and singularity confinement 1 1 7] . This 
study considers ultradiscrete systems defined by a tropically rational map 

(0.1) X n+1 = $(Xn). 

By tropically rational we mean the functions are expressible solely in terms of the binary 
operations ©(max), <8>(+) and 0(— ). We consider X n to be an element of S m , where 
S = KU {- oo}is the max-plus semiring [S], We consider systems of the form (|0. 1 [) to be 
extended cellular automata because the operations of the max-plus semiring allow the of 
the possible restriction of X n to Z m C S m [T5j . This possible restriction of these integrable 
systems to the integers makes them particularly attractive in the sense that, despite their 
simplicity, represent the very core characteristics of integrable systems [TT]. 

The ultradiscretization procedure is a routine procedure that allows one to take a 
discrete dynamical system defined in terms of a single rational mapping, 

(0.2) x n+ i = <f>(x n ), 

which is describable completely in terms of +, X and (subtraction- free), to one of 
the form (|0.ip in a way that is said to preserve integrability[TT]. Ultradiscretization was 
first introduced to relate the Box and Ball svstem[Tl] to integrable lattice equations 16 . 
Various notions of integrability for discrete equations are naturally mapped to analogous 
notions in the semiring setting, such as conserved quantities [TD] and the existence of Lax 
pairs [5]. 

Singularity confinement is widely regarded as a discrete analogue of the Painleve 
property pp. Recently, an analogous notion of singularity confinement has appeared for 
systems of the form (|0.1j) [3l |4]. However, singularity confinement does not naturally 
map via the current form of ultradiscretization to an analogous notion for systems of the 
form([niE]). 
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We present an alternate interpretation of ultradiscretization in terms of a non-archimedean 
valuation field. Recent advances in tropical geometry will allow us to give evidence that 
this new interpretation of ultradiscretization maps singularity confinement from systems 
of the form (|0.2[) to the recently proposed notion of singularity confinement for systems 
of the form (|0.1|) . In §1 we shall present the alternate form of ultradiscretization. In §2 
we shall describe the correspondence between the two forms of singularity confinement. 

1. Extended interpretation of ultradiscretization 

We describe the ultradiscretization procedure as follows: given a subtraction-free ra- 
tional function, f{x), where x = (x%, . . .x rn ) and the S ELI* 6 all positive, we introduce 
corresponding ultradiscrete variables, X = (X\, . . . ,X m ), via the relation Xi — e Xi l t . 
The ultradiscretization of / is the tropically rational function, F, defined by the limit 

(1.1) F(X) = lim elog/Or). 

e->0+ 

To see exactly how this procedure works, let us consider how ultradiscretization applies 
to binary operators: 

(1.2a) X1+X2 — * max(Xi, X2) 

(1.2b) £1x2 - X x +X 2 

(1.2c) xx/x2 -» X 1 -X 2 , 

where — > —00 and 1 — > 0. From the definition given in (11. 1| . it is clear that if one 
considers the ultradiscretization of x\ — x 2 , if X2 > X\, then by the definition above, the 
ultradiscretization requires the logarithm of a negative number. The work of Kasman and 
Lafortune provides methods for overcoming the difficulty presented by the minus sign in 
some cases [6]. 

This correspondence given in (|1.2p motivates the following lemma. 



Lemma 1.1. The ultradiscretization defined by (|l.lj) of any (subtraction-free) rational 
function is equivalent to the replacement of the Xi with Xi and the replacement of binary 
operations defined by (|1.2| . 

Proof. The replacement of the Xi with Xi can be seen from the following calculation 

Xi — > lim eloge Xi / e = lim Xj = Xi- 

To see that the ultradiscretization procedure is equivalent to the replacement of binary 
operations described by (|1.2p . we shall assume the ultradiscretization of f(x) and g(x) is 
F(X) and G(X) respectively. By definition 



f(x)g(x)^ lim elog f(x)g(x) = lim e(\ogf(x)+logg(x))=F(X) + G(X) 

x— >0+ x— »0+ 



and 



f(x)/g(x) -» lim e\ogf(x)/g(x) = lim e(\ogf(x) - logg(x)) = F{X) - G{X) 
x^a+ x— >o+ 

by the properties of the logarithm. This gives correspondences (|1.2b|) and (|1 .2c[l . To see 
(|1.2ap . we note that by definition of the limit, for S > there exists an eo such that 

\elogf(x)-F(X)\<6 
\elogg(x)-G(X)\<6 
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for < e < eo, therefor, under the same assumptions 

e (F(X)-6)/ e K f[x) K e (F(X)+S)/e 
e (G(X)-6)/e < g{x) < e (G(X)+8)/e_ 

By adding these inequalities we have 

e -6/e {e F(X)/e + e G(A)/ e) < f {x) + g{x) < ^(/W/e + e G(X)/«) 

hence 

|elog(/(or) + 3(0;)) - e log(e f + e G(x)/e )| < 5 

for < e < eo- By taking limits as e — * gives the ultradiscrtization of f(x) + g(x) for the 
first term inside the absolute value signs. For the second term inside the absolute value 
signs, since F(X) and G(X) do not depend on e, if F(X) is greater than G(X), then the 
term inside the log is asymptotic to e F ^- x ^ e as e — > + , hence e\og(e F ^ x ^ e + e G ^ x ^ c ) — * 
F(X), and vise versa for G(X). Since (5 > was arbitrary, the ultradiscretization of 
f(x) + g(x) is max(F(X), G(X)). □ 

From this, it is a routine procedure to obtain a mapping of the form (|0.1[) from a 
mapping of the form (|0.2[) . Since our interest lies in integrable cellular automata, we 
consider a simple difference equation, 

(1.3) w n -iw°w n+ i = at n w n + 1 

where |g| ^ 1, t n = toq n and a is constant in n. This equation is considered integrable 
for a G {0,1,2}. Under these conditions, this is a discrete version of the first Painleve 
transcendant^llj. By letting x n — w n -± and y n by w n , we obtain a system of the form 
given by 

Vn 
qt n 

Since the power of x n is an integer, both sides expressible solely in terms of +, x and 
hence is subtraction-free. By replacing a with A, q with Q, t n with T n ,x n by X n and 
y n by Y n , then replacing binary operators on the left and right hand side in accordance 
with (|1.2p . we obtain the ultradiscrete mapping of the form (|0.1[) 

max(A + T n + Y n ,0) - X n - oY n 

T n +Q 

or as a second order equation 

(1.4) W n+1 + aW n + W n -! = max(A + T n + W n , 0) 

where which we now call the ultradiscretized version of the first Painleve transcendant. 

We now present the alternate ultradiscretization as follows. Given a rational function, 
f(x), where x = (x±, . . . , x n ), we let the Xi be algebraic functions in z given by xi — z~ Xi 
where the Xj are rational numbers. Since f{x) is rational in the Xi, f{x) is algebraic in 
z, hence possesses a Puiseux series representation, given by 








(x n 


Yn+1 


-1 


Y n 


Tn+l I 







f(x) = a Q z qa +a x z qi + 
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where qi < qi+i for all i > and ao ^ 0. Given this information, fix) possesses the new 
ultradiscretization, F(X), given by 

(1-5) f{x) - -g 

Lemma 1.2. The ultradiscretization of any (subtraction- free) rational function is equiva- 
lent to the replacement of the Xi with Xi and the replacement of binary operations defined 
by O. 

Proof. Since Xi is rational, Xi is it's own Puiseux series representation, hence, by (jl.5[) . 
Xi — > Xi. We also note that if f(x) — ► and are Puiseux series given by 

f{x) = a z~ F ^ + ... 
g(x) = b Q z- G ^ + ... 

then f(x)g(x) = Z - F (X)-G{X) + _ he nce f( x ) g ( x ) -> F(X) + G{X) and similarly 
f(x)/g(x) -> F(X) - G(X), giving fT2b|) and (fT72c|) . To see (fT2a|) . the addition of the 
two functions is 

fix) + g(x) = a z- F W + bQZ -G(x) + _ _ 
The ultradiscretization of this is - min(-F(X), -G(X)) = max(F(X),G(X)). □ 

This allows us to deduce that the two forms of ultradiscretization agree for subtraction 
free functions. However, in the general case, one is not presented with the same difficulties 
in considering the definition of the ultradiscretization of general (not subtraction-free) 
rational functions. We denote the new ultradiscretization mapping v : C[[t]] — > §, where 
by C [[£]], we mean set of algebraic functions in t. In general, v satisfies the properties 

1. vif) — -co if and only if / = 0. 

2. K/s) = K/) + K<?)- 

3. vif + g)<vif) + vig). 

which implies v is a valuation. However, v also satisfies the following 

3a. vif + g) < max(i/(/), !/(#)). 

This property implies that v is a non-archimedean valuation^- Hence, we have expressed 
the ultradiscretization procedure in terms of a non-archimedean valuation. It is important 
to note that given the two leading terms in the Puissuex series of / and g in / + g do 
not cancel, we have that equality holds in property (3a) of valuations. This allows us 
to map solutions of discrete systems that are not sub-traction free to solutions to their 
corresponding ultradiscrete systems [S]. 

Given a system such as (|1.3[) . one may easily lift the equation to the level of algebraic 
functions in z. We simply make the substitution of w n with W n and constants, {a^}, with 
{z^}. For example, a lifted version of (|1.3p is simply 

^n+l^^n-l - Z A .9 n W n + 1 

where 37~ n +i — 2? n - Given "W n is algebraic in z, by assigning a specific value of z, we 
recover a sequence that satisfies (|1.3[) . Furthermore, given there is no cancellation of 
the dominant powers of z in the calculation of W n so that equality in property (3a) of 
valuations hold for all n, we recover a sequence that satisfies (11.41) . 



^This definition differs siightly from algebraic texts such as Lang[7]. To obtain a non-archimedean 
valuation in the sense of Lang, one simply is required to consider e" rather than v itself. 
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2. Correspondence between the singularity confinements 

Let us first describe the singularity confinement property for ordinary difference equa- 
tions as it is presented originally by Ramani et al.QQ. The singularity confinement prop- 
erty requires that given all singularities, either roots or poles, of a system, that in some 
limiting sense, those singularities do not propagate for all time. We present a simple ex- 
ample to illustrate this definition; we present singularity confinement for the autonomous 
version of equation (|1.3[) . where a = I, to = 1, g = 1 and a = for simplicity. That is, we 
consider the discrete dynamical system defined by the equation 

(2.1) w n -iw n+ i = w n + 1. 

It is clear that we may write this system in the form 




as one can easily see that y n satisfies the same difference equation as w n . Notice that in 
this special case, if wq is free and w± = e, then we encounter a problem in defining u>3 
in the limit as e — > 0. To see whether this system possesses the singularity confinement 
property, we consider the behavior of this system with the above initial conditions. The 
iteration of the initial conditions where wo is free and w\ is e, along with the limit as 
e — > is given by 




The important feature is that despite iterating through a singularity, we are able to retain 
information about the initial condition in the limit as e — > 0. We can sensibly define </> 5 
at (wo, 0) in a limiting sense since the point (wq, 0) is a removable singularity. Hence, we 
say that this singularity is confined. 

Notice that we obtain similar problems in the case where Wq is free and W\ = e — 1 as 
we let e — > 0. We examine the iteration of these initial conditions: 

hm e _ 

(3) - (o 1 ) - (-°0 - (-r.g - - (3) 

The important feature is that (u>o, —1) is a removable singularity of 4> 3 . Another point of 
interest when calculating this map projectively are is the initial conditions (wq, oo) which 
can be seen by letting e — > oo in the first example. In all these examples, it should be 
clear that the singularities are confined, since </> 5 is equivalent to the identity map. 

Definition 2.1. A mapping, (j> : S — > S, is said to have the singularity confinement 
property if for every singularity, p € S, there exists an n such that p is a removable 
singularity of <f) n . 
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One key point in the above example is that (j> 5 is equivalent to an isomorphism, despite 
the fact that (j) 3 is not. This is seen by the way that <fi 3 sends (wq, 0) to (oo, oo) regardless 
of the value of wq- We may say that a whole copy of Pi has been sent to (oo, oo). 

In this case, we may consider the blow-up of Pf at points Pi = (— 1,0), P2 = (0, — 1) 
and P3 = (00, 00) resulting in the surface ScPj, given by 

J x+1 x y 

6 = i {x,y,Pi = ,P2 = — -~r,P3 = ~ 

I y y+ 1 x 

If we lift <p to the mapping <f>' : S — ► S, (j)' is an isomorphism. The resulting surface 
is a highly degenerate (and simplified) form of the much celebrated surface of initial 
conditions [15]. 

We describe the ultradiscrete analog of singularity confinement in a similar manner to 
discrete mappings. The ultradiscrete singularity confinement requires that any points of 
non-differentiability do not propagate for all time. We present the ultradiscrete of the 
same simple example given above. We consider the mapping defined by the iterative 
scheme given by 

(2.2) W n -x + W n +i = max(W n , 0). 

This may also be viewed as the 2 dimensional system defined by the mapping 

(2 ' 3) {r n ti ) = * ( W ) = (max(0, Y n ) - X n 

The mapping $ possesses a point on non-differentiability at the point corresponding 
to initial conditions where Wq is free and Wi = 0. Notice that although W n is not 



n 


W n : S = 


W n : 5 > 


S < 





W 


Wo 


Wq 


1 


5 


5 


5 


2 


-Wo 


5 -Wo 


-Wo 


3 


max(0, -Wo) 


max(— 8, — Wo) 


max(0, -Wo) - S 


4 


max(0, Wo) 


max(0, Wo) - S 


max(0, Wo) - <5 


5 


Wo 


Wo 


Wo 


6 


5 


5 


6 



Table 1 . We compute the iterates of W„ when Wo is free and W n = S 
where \S\ is a assumed to be small. This allows us to calculate the left 
and right derivatives at W\ =0. 



differentiable at n = 2 and 3 at Wi = 0, W4 is differentiable at W\ = 0. That is that 
points of non-differentiability do not propagate for all time. 

Definition 2.2. A mapping, $ : S — > S, is said to have the ultradiscrete singularity 
confinement property if, for every point, p, such that $ is not differentiable at p, there 
exists an n such that <I>™ is differentiable at p. 

Just as in the discrete case, the singularities of (|2.3[) arc confined in the ultradiscrete 
sense since $ 5 is the identity map. 
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Lemma 2.3. Given a subtraction free rational function, f(x) = , with ultradis- 
cretization, F{X) = Fq(X) — F\{X), the image of the roots of fi and fi are points of 
non- differentiability of Fq and Fi. 

Proof. We simply remark that this is a reformulation of Theorem 2.7 of [12]. Given a 
subtraction free polynomial, fo, the ultradiscretization, Fq, the algebraic function in t 
provided in the proof of Theorem 2.7 of [12] is the algebraic function in z considered in 
the alternate ultradiscretization obtained by replacing the variables at and x with z Ai 
and z x respectively. It is clear that the points of non-differentiability of the valuation 
considered in |12j is in one-to-one correspondence with the points of non-differentiability 
of the valuation considered in the alternate ultradiscretization in light of the relation 
max(Xi, . . . , X n ) = — min(— X\, . . . , —X n ). Finally, by the properties of the valuation, 

v{f{z x )) = Hfo(zl) - v{f x {z x )) = F - Fx, 

hence the points of non-differentiability is the union of the points of non-differentiability 
of .Fo and Fi. □ 

We examine the effect of the ultradiscretization of singularity confinement of over the 
system defined by (|2.1|) after lifting the system to the field of algebraic functions in z and 
applying valuation. That is we consider the system defined by 

(2.4) #n+l#n-l ~ W n + 1- 

where the variables W n are rational functions of z. To consider singularity confinement, 
we have a few points of interest: Wx — z~°° — and Wx = — z° = — 1 where in both 
cases, we will consider Wo to be of the form z w ° where Wo is free. 



n 


w n 


W n 


Root/Poles in Wx 


ND Points in Wx 





w 


W 










1 


Wi 


Wi 




{0} 


{-co} 


2 


% 


max(0,V7i) 


-Wo 


{-1} 


{0} 


3 
4 


1 + Wo+Wi 


max(0,T^ , W x ) - 
max(0, Wo) 


-Wq-Wx 
- Wx 


{-1-^0,0} 
{0} 


{max(0, W ), -oo} 
{-oo} 


5 




W 










6 




Wi 




{0} 


{-co} 



Table 2 . The iterates of W n and W n as well as the poles and roots of W n 
as a function of W\ and the points of non-differentiability (ND Points) 
as a function of Wx ■ 



By considering table [21 one is use Lemma 12.31 to determine when one would expect 
points of non-differentiability from the Roots and Poles of W n as functions of Wx alone. For 
example, at Wx = 0, we see this is v(—l), hence one expects that W2 is not differentiable 
at as a function of Wx- Also if Wq < 0, then W3 will also not be differentiable in W\ 
since if Wo < then v(l + Wo) = 0. In the case Wo = —00, we consider this to be ^(0), 
hence we expect this to not be integrable at n = 1, 3 and 4. 

In the first case, we consider an initial condition where Wx is close to 0, namely an 
initial condition where Wx — z~ L where we consider Wx — > as L — * 0. Hence, we may 
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consider the valuation under the assumption that L 0. 



n 


3T„ 


W n = u{W n ) : L » 







W 


1 


Z" L 


-L 


2 


z -Wo(l + z -I>) 


-W 


3 


z- w °(l + z L + z L+w °) 


L - W + max(Wo,0) 


4 


z L (l + z Wo ) 


L + max(0,W ) 


5 


z w 


W> 


6 


z~ L 


-L 



The first feature of this is that due to the subtraction-free nature of W n as rational 
functions of z with these initial conditions W n — v(W n ) satisfies (|0.1j) . The second 
feature is that, as predicted by Lemma l2~3l the mapping is not differentiable at n = 1, 3 
and 4 as the functions diverge. Also notice that by assigning a value to z, as L — > oo, 
z~ L — > or oo. In either case, what is observed is singularity confinement for (|0.2[) . 

We now turn to the initial condition around W\ — —z°, given by — z where \8\ <C |Wo|. 
There are two cases we are required to consider: the case in which 8 < and the case 
where 8 > 0. The valuation of W n for each n depends on 8. 



n 




W n 


= v(W n ) : 8 > 




= v(Wn) ■ 8 < 





Z W 




Wo 




Wo 


1 


-z s 




8 




8 


2 


z- Wo (l-z 5 ) 




-Wo 




-8 -Wo 


3 


z- s - w °(z s -l-z w °) 


-Wo 


+ 8 + max(0, W ) 


-Wo 


+ max(0, Wo + 8) 


4 


-z- s (l + z Wo ) 


8 


+ max(0, Wo) 


8 


+ max(0, Wo) 


5 


Z W 




Wo 




Wo 


6 


-z 5 




8 




8 



The first feature, as was observed in the previous case, is that W„ = v{W n ) satisfies (|0.ip 
for cases in which 8 > and 8 < 0. This can be guaranteed by sufficiently generic (and 
generally sufficiently small) values of 8 so that there is no cancellation of highest power 
terms of z in the calculation of W n . Furthermore, as predicted by the use of Lemma [2751 
the function is not differentiable at n = 2, while for n = 3, W n is differentiable if Wo < 
and is not differentiable if Wo > 0. 

As remarked at the end of §2, the correspondence between the two types of singular- 
ity confinement can be summarized as follows. After lifting the equation to the level of 
algebraic functions in a variable z, by assigning a specific value of z, we recover singu- 
larity confinement as it is defined for discrete mappings. However, given one takes initial 
conditions sufficiently and generically close to the singularities, one recovers ultradiscrete 
singularity confinement via the application of the valuation. 

3. Discussion 

We presented a very simple autonomous example to demonstrate the correspondence 
between the two types of singularity confinement. We could have chosen any one of the 
equations presented in [TT], however, the system chosen was one of the simplest systems 
we could find that displayed all the desired characteristics. It is worth noting that al- 
though this correspondence works for this example, and indeed for all the ultradiscrete 
Painleve equations that featured in [TT], what we have demonstrated is not a one-to-one 
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correspondence between the two singularities. In fact the inverse of the valuations of 
hopelessly multi- valued. A case in which we have not ruled out is one in which by iter- 
ating out of one singularity, you are in the image of another pole or root. Furthermore, 
we have provided a surface of initial conditions in which the mapping <j> was an isomor- 
phism, however, further study is required to determine what the analogous notion is for 
ultradiscrete systems such as $. 
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